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Complex- Temperature Singularities of Ising Models 
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We report new results on complex-temperature properties of Ising models. These include studies of the s = 1/2 
model on triangular, honeycomb, kagome, 3T2^, and 4-8^ lattices. We elucidate the complex-T phase diagrams of 
the higher-spin 2D Ising models, using calculations of partition function zeros. Finally, we investigate the 2D Ising 
model in an external magnetic field, mapping the complex-T phase diagram and exploring various singularities 
therein. For the case PH = j7r/2, we give exact results on the phase diagram and obtain susceptibility exponents 
7' at various singularities from low-temperature series analyses. 



1. Introduction 

During the past year, in collaboration with Vic- 
tor Matveev, we have carried out a number of 
studies of complex-temperature (CT) properties 
of Ising models We report on some results 

from these here (see also There are several 

reasons for studying the properties of statistical 
mechanical models with the temperature variable 
generalized to take on complex values. First, one 
can understand more deeply the behavior of vari- 
ous thermodynamic quantities by seeing how they 
behave as analytic functions of complex temper- 
ature. Indeed, complex-T singularities have im- 
portant influences on the behavior of such func- 
tions on the real K axis (see definition in eq. (H)). 
Second, one can see how the physical phases of a 
given model generalize to regions in appropriate 
complex-temperature variables. Third, a knowl- 
edge of the complex-temperature singularities of 
quantities which have not been calculated exactly 
helps in the search for exact, closed- form expres- 
sions for these quantities. This applies, in par- 
ticular, to the susceptibility of the 2D zero- field 
Ising model, and to the Ising model for higher 
spin or external magnetic field, neither of which 
has been solved above d = 1. Some of the ear- 
liest works on complex-T properties of the Ising 
model include Refs. For comparison, a 

different direction is to study the case in which 
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the temperature is kept physical but the exter- 
nal field H is complex We also report new 
results which both K and H are complex. 

2. Model 

The Ising model, and our notation, are spec- 
ified as follows: the partition fimction is Z = 

Eia^yP'"^ with 



(1) 



where at = ±1, /3 = {ksT)"^ , and H denotes a 
possible external field (zero unless otherwise in- 
dicated). We use the standard notation 



K = (3J , h = (3H 



/i = e 



-2h 



(2) 
(3) 



V = tanhii', and the elliptic modulus variable 
fc< = 4u/(l — uY . We define the reduced free 
energy as f — —/3F = limTv^oo \nZ and the 
reduced susceptibility as x = P~^X- 

3. Phase Diagrams (Homopolygonal Lat- 
tices) 

Owing to the symmetries under the shifts K — > 
K + niTT and K K +{2n+l){iTT/2) g, there is 
an infinite replication of complex phases as func- 
tions of K; this is reduced to a single set by using 
the variables v or z. For a bipartite lattice, the 
phase boundaries are invariant under z — > — z. 
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so that the most compact way to show them is 
in terms of the variable u = z 2. In Ref. 
we showed that for the square lattice the phase 
boundaries in the u plane are given by a Pascal 
limagon. The inner and outer branches of this 
limagon cross at the point u = —1. 

We have extended this analysis to the triangu- 
lar and honeycomb lattices |^]. As in the case of 
the square lattice, the continuous locus of points 
where / is non-analytic can be determined ex- 
actly. (/ is also trivially non-analytic at the 
points \K\ — oo, but since these do not sepa- 
rate any phases, they will not be important here.) 
In contrast to the square lattice, this (continu- 
ous) locus of points includes arcs or line segments 
which protrude into, and terminate within, cer- 
tain phases and hence do not completely separate 
any phases. For example, for the triangular lat- 
tice, in the u plane, the above locus of points con- 
sists of the union of the circle u = (—1-1- 2e'^)/3 
and the semi- infinite line segment —oo<u< 
— 1/3. The interior (exterior) of the above cir- 
cle forms the FM (PM) phase. M diverges at 
u — —1/3 with exponent (3 = —1/8 and is discon- 
tinuous across the line segment —1 < u < —1/3. 
We have proved that such a divergence in M au- 
tomatically implies a divergence in x Q . 

4. Heteropolygonal Lattices 

While the square, triangular, and honeycomb 
lattices each involve only one type of regular poly- 
gon (and thus may be denoted "homopolygonal" ) , 
there are also regular tilings of the plane (het- 
eropolygonal lattices) which involve more than 
one type of regular polygons, but still have 
all vertices equivalent, and all bonds of equal 
length. We have determined the exact complex- 
temperature phase diagrams of the 2D Ising 
model on three regular heteropolygonal lattices: 
3 • 6 • 3 • 6 (= kagome), 3 • 12^, and 4 • 8^ (bath- 
room tile), where the notation denotes the regu- 
lar n-sided polygons adjacent to each vertex 
We have also worked out the exact complex- 
T singularities of M. In particular, we found 
the first case (the 4 • 8^ lattice) where, even 
for isotropic spin-spin exchange couplings, the 



nontrivial non-analyticities of the free energy of 
the Ising model lie in a two-dimensional, rather 
than one-dimensional, algebraic variety in the 
z = e~^^ plane. 

5. Complex T Singularities 

Complex-T singularities exhibit a number of 
new features not observed for physical critical 
points. First, as we showed by exact results 0, 
universality, in the sense of lattice-independence 
of critical exponents, is violated at u = — 1, 
since the magnetization exponent /3 is 1/4 for the 
square lattice but 3/8 for the triangular and hon- 
eycomb lattices. Second, on the square lattice, 
we proved ||] that as one approaches the bound- 
ary of the complex-T extension of the PM phase, 
X has at most finite non-analyticities except for 
its divergence at the physical critical point Uc, 
and hence, in particular, at m = — 1 = u^, it fol- 
lows that 7s < 0. Analyses of low-temperature 
series indicate that x diverges at u = — 1 as this 
point is approached from within the FM phase 
[|l|,|l| (and AFM phase [|l|), with the exponent 
7s numerically consistent with 3/2. This value 
also follows from the relation 7^ = 2(7' — 1) 
Hence, 7s 7s, which is unprecedented for phys- 
ical critical points. Third, the correlation length 
exponent v'^ has different values depending on 
whether one extracts it from non-diagonal cor- 
relation functions or diagonal correlation func- 
tions. Several scaling relations are violated; for 
example, + 2/3s -I- 7s ^ 2 for the singularity 
at M = — 1 as approached from within the PM 
phase. For honeycomb lattice, we have shown |^ 
that a' + 2(3 + j'^2 (the RHS is consistent with 
being equal to 4) at the CT singularity at z = — 1, 
as approached from within the FM phase, and 
have discussed the origin of this. 

6. Higher-Spin Ising Model 

We have carried out a study of the complex-T 
properties of higher-spin 2D Ising models Q . The 
method which we use for this purpose is the cal- 
culation of complex-T zeros of the partition func- 
tion on finite lattices. It is known from exactly 
solved cases that as the lattice size increases, the 
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CT zeros of Z occur on, or progressively closer 
to, the continuous locus of points where / is non- 
analytic in the thermodynamic limit. Hence, by 
studying these zeros, one can make reasonable in- 
ferences about this locus of points and the cor- 
responding complex-T phase diagram. Specific 
results are presented for s = 1, 3/2, and 2 in 
Ref. 1^. We infer several general features of the 
phase diagrams and motivate a conjecture that 
the number of arcs protruding into the complex- 
temperature extension of the FM phase is 4[s^] — 2 
for s > 1, where [a;] denotes the integral part 
of X. From comparison with an analysis of low- 
temperature series expansions further sup- 
port is obtained for our conjecture that the mag- 
netization diverges at endpoints of such arcs ter- 
minating in the FM phase. We also present an 
exact determination of the complex-temperature 
phase diagram for the ID spin s Ising model. We 
find that in the complex Us = e^^/'* plane (i) 
it consists of A^c.iD = 4s^ infinite regions sep- 
arated by an equal number of boundary curves 
where the free energy is non-analytic; (ii) these 
curves extend from the origin to complex infinity, 
and in both limits are oriented along the angles 
On = (1 + 2n)7r/(4s2), for n = 0, ...,45^ - 1; (iii) 
there is a boundary curve (line) along the nega- 
tive real Us axis if and only if s is half-integral. 
An interesting connection is noted between some 
features of the ID and 2D complex-temperature 
phase diagrams for higher-spin Ising models. 



7. Finite Magnetic Field 

Although the Ising model has never been solved 
in an arbitrary external magnetic field, Lee and 
Yang [Q did succeed in solving exactly for the 
free energy and magnetization for the model on 
the square lattice for f3H = h = i-K/2. This solu- 
tion yields further insight into the structure of the 
model. It is interesting to investigate the proper- 
ties of this model when one also generalizes the 
temperature to complex values, and we have done 
this . We have first determined the continuous 
locus of points where / is non-analytic for this 
model: in the u plane this is the union of (i) the 
unit circle and (ii) the finite line segment along 



the negative real axis 
1/ue < u < Ue 
where 

Zie = -(3-2-3/2) 



(4) 



(5) 



The complex-T phase diagram of this model con- 
sists of a FM phase inside of, and an AFM phase 
outside of, the unit circle in the u plane. We find 
that the specific heat C has the following singu- 
larities: divergences, with exponent a'g = 1 at Me 
and 1/ue, and finite logarithmic non-analyticities 
sA u = —1 = Us and u = I (a'j, = 0, a'l = 0). The 
uniform and staggered magnetizations M and 
Mst vanish with exponent /3s = 1/2 at u = — 1 
and diverge at it = 1 with exponent j3i — —1/4. 
M and M^t diverge, respectively, at and 1/we, 
with exponent /3e — —1/8. 

From the Baxter-Enting low-temperature, 
high-field series expansion of the partition func- 
tion, we have extracted the low-temperature se- 
ries for the susceptibility x to 0{v?^). Analyzing 
this series, we find that x has divergent singulari- 
ties (i) at u = Ue with exponent 7^ = 1.25 ±0.01; 
(ii) at w = 1, with exponent 7^ — 2.50 ± 0.01; 
and (iii) at it = Us = —1, with exponent 7^ = 
1.00 ± 0.08. We infer the exact values of these 
exponents to be 
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(6) 

(7) 
(8) 



We have subsequently determined the exact 
complex-T phase diagram for the model aX, h = 
z7r/2 on the honeycomb {he) and triangular [t) 
lattices 0. For the t lattice, we find that this con- 
sists of a single (FM) phase, and the continuous 
locus of points where / is non-analytic consists of 
the union of (i) the circular arc 



(9) 



where 9ce — arctan{4:^/2/7) ~ 38.9°, correspond- 
ing to the endpoints Uce = e 



and where 
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Uce = (~1 + 2'^/^i)/9, and (ii) the semi- infinite 
line segment 



oo < It < -1/2 



(10) 



For the he lattice, we have proved that all of our 
results on the complex-T phase diagram and the 
singularities of x for /i = on the he lattice also 
apply for /i = i7r/2 with the replacement of z by 

— z. 

In Ref. we have recently investigated the 
complex-T phase diagram of the square-lattice 
Ising model for nonzero external magnetic field 
H in the range —oo < H < oo, i.e. < fj, < oo 
where /i is given by (^. We have also carried 
out a similar analysis for ~oo < n < 0. Without 
loss of generality, one may restrict to the range 

— 1 ^ M ^ 1- The results for this interval provide 
a new way of continuously connecting the two 
known exact solutions of this model, viz., that at 
H = 1 {h = 0) and fj, = —1 [h — i7r/2). Our meth- 
ods include calculations of complex-temperature 
zeros of the partition function and analysis of low- 
temperature series expansions. For real nonzero 
i?, we find that the inner branch of a limagon 
bounding the FM phase breaks and forms two 
complex-conjugate arcs which move away from 
the real u axis as increases. We have studied 
the singularities and associated exponents of ther- 
modynamic functions at the endpoints of these 
arcs. For < 0, there are two line segments of 
singularities on the negative and positive u axis, 
and we carry out a similar study of the behavior 
at the inner endpoints of these arcs, which consti- 
tute the nearest singularities to the origin in this 
case. 
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8. Conclusions 

Thus, the study of complex-temperature singu- 
larities of the Ising model provides many interest- 
ing insights. A number of topics deserve further 
investigation, such as related spin models and be- 
havior in d > 2. 
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